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摘 要 : 冰川 覆盖 流域 的 雪 冰 融 水 对 河川 径流 有 重要 的 调节 作用 。 气 候 变 化 影响 雪 冰 融 水 过 程 和 
数量 变化 ,河川 径流 过 程 和 径流 量 相应 变化 ,其 程度 与 流域 冰川 情况 相关 。 通 过 利用 CMIP5 气候 
模式 输出 气象 数据 驱动 流域 水 文 模型 ,模拟 研究 天 山地 区 3 个 不 同 冰 川 覆盖 率 河 流 ( 库 玛 拉克 河 、 
玛 纳 斯 河 、 库 车 河 ) 的 径流 对 气候 变化 的 响应 。 结 果 表 明 : 随 着 未 来 气温 和 降水 的 持续 增加 ,3 个 流 


域 的 雪 融 水 均 有 增加 , 冰 融 水 变化 受 冰川 覆盖 鱼 


积 的 影响 ,在 各 个 流域 变化 不 一 致 。 径 流 变化 主 


要 受降 水 增加 和 雪 冰 融 水 变化 的 综合 影响 ,在 未 来 情景 下 各 流域 径流 均 有 增加 ,分 别 增加 了 35.8%-~ 
14.3%( 库 车 河 )、.2.9%~11.49%( 玛 纳 斯 河 )、12.9%~47.19%( 库 玛 拉 克 河 ), 且 冰川 覆盖 率 越 大 的 流域 ， 
预 估 径流 不 确定 性 变化 区 间 受 冰 融 水 影响 越 大 。 预 估 3 个 流域 的 径流 、 雪 冰 融 水 年 内 分 布 变化 表 
明 , 各 河流 的 春季 融雪 时 间 提 前 和 融雪 量 增加 使 得 流域 春季 径流 量 较 历 史 时 期 增 大 ;在 夏季 , 受 雪 


冰 融 水 变化 的 影响 库 车 河 、 玛 纳 斯 河 夏季 径流 峰值 


量 减 小 ,而 库 玛 拉克 河 径流 峰值 量 增 加 , 且 预 估 


的 各 流域 夏季 径流 变化 不 确定 性 区 间 明 显 大 于 其 它 季 节 。 
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IPCC 第 五 次 气候 评估 报告 (AR5 ) 指 出 , 近 20 a 
来 全 球 山 地 冰川 整体 处 于 退缩 状态 ,北半球 的 积 雪 
面积 持续 缩小 "一 , 且 有 冰川 覆盖 的 流域 其 径流 对 气 
修 变 化 响应 更 为 敏感 时。 天山 广泛 分 布 冰川 ,是 众 
多 流入 中 国 新 疆 和 中 亚 五 国 河流 的 发 源 地 "1, 是 下 
游 地 区 的 人 类 生存 和 发 展 的 主要 水 源 "。 一些 研究 
者 "分 析 了 天 山 观 测 冰川 、 各 水 系 或 流域 冰川 过 
去 几 十 年 间 的 变化 ,表明 受气 候 变 暖 影响 天 山 不 同 
区 域 已 发 生 不 同 程度 的 退缩 。 冰 川 融 水 虽然 可 以 
在 短期 内 提供 丰富 的 水 量 , 但 是 其 总 储量 在 减少 ， 
如 果 冰 川 消融 继续 加 速 , 未 来 以 融 冰 和 融雪 为 主要 
MAS B vb RT Bez HH BUR BK vH UL. BG 
着 气候 持续 变化 , 越 来 越 多 的 学 者 利用 不 同方 式 模 
拟 的 未 来 气候 情景 数据 驱动 水 文 模型 模拟 预 估 径 
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流 变化 "1。 

有 具有 冰川 模块 的 水 文 模型 是 研究 雪 冰 流域 水 
文 过 程 变 化 特征 的 一 种 有 效 工具 ,通过 模型 模拟 可 
以 很 好 的 表现 出 冰川 消退 与 流域 产 流 的 关系 "" ,其 
中 气候 模式 输出 的 气候 变化 情景 数据 驱动 水 文 模 
型 是 预 估 未 来 气候 变化 对 水 文 过 程 影响 的 常用 方 
法 。 近 年 来 已 有 较 多 研究 利用 CMIP5(Coupled 
model intercomparison project phase 5 ) 气 候 模 式 数 据 
驱动 水 文 模型 模拟 预 估 中 亚 天 山 、 咯 喇 昆 仑 山 、 青 
藏 高 原 .喜马拉雅 山上 的 一 些 雪 冰 流 域 未 来 几 十 年 
径流 变化 趋势 . 积 雪 融 水 变化 及 冰川 面积 变化 ””。 
Bosshard 55 ^* I Vetter 等 ”分别 采用 5 个 和 10 个 气 
候 模式 预 佑 Rhine 流 域 未 来 径流 变化 ,得 到 了 相同 
的 径流 变化 趋势 但 变化 幅度 存在 差异 。Zhang 等 " 
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选取 了 5 个 气候 模式 驱动 水 文 模型 模拟 预 估 了 乌 鲁 
木 齐 河 未 来 几 十 年 径流 的 变化 ,结果 显示 冰川 融 水 
在 未 来 会 呈 先 增加 后 减 小 的 趋势 ,气候 模式 输入 的 
气象 数据 对 径流 模拟 影响 较 大 。CMIP5 发 布 了 50 
多 个 气候 模式 , 且 选 取 气 候 模 式 越 多 , 预 佑 径流 变 
化 的 不 确定 性 区 间 越 大 ”。 本 文采 用 CMIP5 中 能 
够 获得 的 39 个 气候 模式 ,这 些 气 候 模式 模拟 的 气温 
和 降水 输入 水 文 模型 可 以 给 出 流域 水 文 过 程 模拟 
预 佑 结果 的 变化 范围 。 

在 雪 冰 流域 , 积 雪 和 冰川 作为 降水 最 终 转化 为 
径流 过 程 中 的 一 种 中 间 状 态 ,对 降水 .气温 变化 的 
响应 敏感 ™'。 一些 研 究 结果 表明 ,过 去 几 十 年 间 天 
山 各 地 区 径流 变化 的 主要 影响 因素 是 不 同 的 ”。 
例如 ,高 前 兆 等 ”对 比 了 阿克苏 河和 开 都 - 孔 霍 河 
气温 、 降 水 和 径流 的 关系 ,结果 表明 开 都 河 的 径流 
变化 幅度 要 大 于 阿克苏 河 ,并 指出 以 冰川 融 水 补给 
为 主 的 河流 比 以 降水 补给 为 主 的 更 加 稳定 。 未 来 
几 十 年 ,不 同 冰川 覆盖 率 的 流域 对 气候 变化 的 响应 
又 是 如 何 呢 ? 基于 以 上 问题 ,利用 39 个 气候 模式 驱 
动 具 有 模拟 冰川 水 文 过 程 能 力 的 水 文 模型 SWAT- 
RSG ” ,探讨 天 山地 区 冰川 覆盖 率 不 同 的 3 个 流域 


(a) 研究 区 冰川 覆盖 图 


(c) 库 车 河流 域 
0 2500km 
i 


高 程 /m 


协和 拉 水 文 站 Fm 


1289 兰 王 水 文 站 


其 径流 对 气候 变化 的 响应 ,同时 也 可 以 为 该 地 区 合 
理 优 化 水 资源 配置 提供 相应 的 理论 依据 。 


1 研究 区 数据 与 方法 


1.1 研究 区 概况 

天 山 山脉 (39°~45°N , 699-95? E) Fe VG BR ES npn Fs] 
新 疆 、 哈 萨克斯 坦 吉尔吉斯 斯 坦 和 乌 效 别克 斯 坦 4 
国 ,整个 天 山 可 被 划分 为 中 天 山 西天 山东 天 山 和 
北 天 山 ""5。 本 文选 取 天 山地 区 冰川 覆盖 率 不 同 的 
库 玛 拉克 河 ` 玛 纳 斯 河和 库 车 河 的 产 流 区 作为 代表 
性 流域 ,人 研究 气候 变化 对 雪 冰 水 文 过 程 .径流 过 程 
4 影响 ;流域 位 置 和 冰川 覆盖 情况 见 图 1, 特征 见 表 
1。 库 玛 拉 克 河 流域 冰川 面积 覆盖 率 最 大 为 21.0%， 
径流 主要 由 雪 冰 融 水 补给 为 主 。 玛 纳 斯 河流 域 冰 
川 面积 覆盖 率 约 为 10.0% ,径流 主要 由 夏季 降雨 和 
雪 冰 融 水 补给 为 主 。 库 车 河 冰 川 面积 最 小 ,冰川 覆 
xt AUR 0.596 ,其 径流 补给 主要 来 源 于 夏季 降雨 。 
1.2 数据 

SWAT-RSG 模 型 构建 所 需要 的 基础 数据 包括 : 
流域 数字 高 程 模型 (DEM) .土壤 、 土 地 覆盖 、 冰 川 、 
气象 和 径流 数据 等 。 


(d) 玛 纳 斯 河流 域 
0 2500 km 


4923 
= 855 


图 1 库 玛 拉克 河 、 玛 纳 斯 河和 库 车 河流 域 位 置 及 冰川 覆盖 示意 图 


Fig. 1 Location of Kumaric River Basin, Manas River Basin, Kuqa River Basin and the distribution of glacier 
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表 1 研究 区 流域 面积 冰川 面 积 、 高 程 和 气候 特征 统计 


Tab.1 Statistics of basin area, glacier area, elevation, and climatic features in study area 


河流 水 文 站 。 流域 面积 hua — — 情况 EE - a 
fi f/km? LEE A 平均 最 低 最 高 气温 /°C 降水 /mm 
库 玛 拉克 河 办 和 拉 12990 2725 21.0 3732 1434 7070 10.8 130.0 
马 纳 斯 河 肯 斯 瓦特 5163 493 9.5 3252 846 5145 7.2 215.0 
库 车 河 兰 干 2912 15 0.5 2515 1271 4488 11.4 66.5 


n 


DEM 数据 采用 ASTER GDEM (ASTER global 
digital elevation model, http://gdem.ersdac.jspacesys- 
tems.or.jp/index.jsp) ,分辨 率 为 90 mo TE ArcSWAT 
中 利用 DEM 数据 生成 子 流域 边界 、 河 道 、 高 程 带 和 
坡度 等 。 土 壤 数 据 从 HWSD (Harmonized world soil 
database ) 数 据 库 提取 。 提 取 数 据 包括 土壤 名 称 分 
层 数 质地、 容重 等 。 而 土壤 有 效 水 含量 、 饱 和 水 力 
传导 度 等 则 需要 通过 SPAW (Soil plant atmosphere 
water) 软 件 进行 处 理 计算 。 根 据 以 上 方法 得 出 的 土 
壤 参 数 写 和 人 到 SWAT 模 型 中 的 用 户 土壤 数据 库 中 。 
土地 利用 图 采用 欧洲 航天 局 2009 年 新 版 世界 卫星 
地 图 GlobCover2009 ,空间 分 辩 率 为 300 m, 

冰川 数据 主要 来 源 于 中 国 第 一 次 冰川 编目 数 
据 集 (http:/westdc.westgis.ac.cn ) ,第 二 次 冰川 编目 数 
据 集 ' 半 。 库 玛 拉克 河流 域 产 流 区 大 部 分 位 于 吉尔 
吉 斯 斯 坦 , 冰 川 数 据 采 用 AsiaCryoWeb 发 布 的 TS- 
GV1(Tien Shan Glaciers V1.0)55 数据 集 。 冰 川 信 息 
包括 地 理 坐 标 面积. 长度、 延伸 方向 、 始 末端 高 程 
和 平均 高 程 以 及 冰川 类 型 和 冰 磺 类 型 ,冰川 矢量 边 
界 数 据 。 以 上 冰川 信息 作为 GHRU 方 法 建 模 的 输 
入 数据 。 径 流 数据 来 源 于 《中 国 水 文 年 鉴 ) 统 计 的 
月 值 数据 , 玛 纳 斯 河 的 肯 斯 瓦特 水 文 站 为 1961— 
2005 年 ; 库 车 河 的 兰 干 水 文 站 为 1961 一 1990 年 ; 库 
玛 拉 克 河 的 协和 拉 水 文 站 为 1961 一 2005 年 。 

3 个 流域 的 气象 观测 数据 情况 。 由 于 气象 观测 
站 通常 位 于 海拔 较 低 的 位 置 ,同时 流域 内 气象 观测 
站 点 稀少 ,高 海拔 地 区 缺乏 观测 数据 ,研究 中 进 一 
步 采 用 网 格 气象 数据 APHRODITE (Asian precipitation 
highly-resolved observational data integration towards eva- 
luation of water resources ) 5 的 降水 数据 和 PCMFD 
(Princeton' s global meteorological forcing dataset ) ^" 
的 气温 数据 。 这 2 个 数据 集 已 在 中 亚 天 山 wl EG 
仑 山 .青藏 高 原 等 缺乏 长 期 观测 气候 资料 的 高 海拔 
地 区 得 到 很 好 的 应 用 ”。 


注 : 冰 川 面 积 统计 来 自 于 中 国 第 二 次 冰川 编目 数据 集 , 时 间 2007s2 光 和 AsiaCryoWeb 发 布 的 Tien Shan Glaciers V1(TSGV 1) ,时 间 2002s55 。 


1.3 水 文 模型 

SWAT-RSG 模型 适用 于 模拟 冰川 ,流域 .区域 多 
尺度 上 的 气候 变化 对 雪 冰 水 文 过 程 和 径流 的 影响 ， 
在 天 山 WAAR AR . 喀 喇 昆仑 山 等 雪 冰 流域 的 应 用 中 
取得 了 较 好 效果 ,本文 将 利用 该 水 文 模型 模拟 
研究 区 流域 水 文 过 程 。 该 模型 是 Luo 等 ”将 冰川 水 
X Me] JW AA FG (Glacier hydrologic response units, 
GHRUs ) fil A $ SW AT [5207 B , ft FEL A DL UN. 
度 融 雪 、 融 冰 和 降雨 产 流 过 程 的 模拟 能 力 ,进而 对 
雪 冰 流域 径流 模拟 。GHRUs 方 法 ,计算 单条 冰川 物 
质 平衡 ,并 用 冰川 体积 -面积 尺度 关系 (Volume-area 
scaling relation ) 模 拟 冰 川 动态 变化 ;冰川 物质 平衡 
计算 中 ,积累 考虑 了 雪 转 冰 、 融 水 再 冻结 ,消融 考虑 
了 升华 ` 冰 上 和 雪 融 水 、 冰 融 水 ; 雪 、 冰 融 水 采用 度 日 
因子 算法 ;冰川 径流 由 冰 上 和 雪 融 水 、 冰 融 水 `. 冰 上 两 
产 流 3 部 分 构成 ;每 个 GHRU 划分 高 程 带 , 以 适应 降 
水 .气温 随 高 程 的 变化 。 方 法 的 详细 论述 见 文献 
Luo SP". 
1.4 模型 模拟 效果 评价 指标 

模型 率 定 和 验证 阶段 ,模拟 的 径流 量 和 实测 值 
的 比较 采用 NSE (Nash-sutcliffe efficiency ) ^^ fll PBIAS 
(Percent bias) 52 4-48 bs YE ffr; NSE 主要 评价 模拟 
的 径流 过 程 和 实测 值 的 吻合 情况 ;PBIAS 主要 评价 
模拟 和 实测 径流 量 的 数量 偏差 。 采 用 Moriasi 5 
给 出 划分 模型 参数 化 优 劣 的 标准 进行 评价 , 见 表 
2。NSE 和 PBIAS 计 算 方 法 如 下 : 


NSE = 1 - | -——— ———— (1) 


PBIAS= 1 - — (2) 


式 中 : 0” 为 第 ;次 的 观测 径流 量 ; O” 为 第 ;次 的 模 
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X2 NSE 和 PBIAS 对 径流 模拟 能 力 评 价 
Tab.2 Evaluation of simulated streamflow by nash- 
sutcliffe efficiency (NSE) and percent bias (PBIAS) 


NSE PBIAS/% 结果 评价 
NSE<0.50 PBIAS>+25 不 满意 
0.50<NSE<0.65 +15<PBIAS<+25 满意 
0.65<NSE<0.75 +10<PBIAS<+15 好 
0.75«NSEx1 PBIAS<+10 非常 好 


拟 径流 量 ; 0 为 观测 径流 量 的 均值 ;为 径流 序列 
1.5 气候 模式 选取 与 统计 降 尺 度 

全 球 气候 模式 是 研究 气候 变化 机 制 和 预 估 未 
来 气候 变化 对 水 文 过 程 影响 的 重要 工具 “。CMIP5 
给 出 了 20 多 个 模式 研究 中 心 发 布 的 50 多 个 气候 模 
式 在 RCP2.6、RCP4.5 .RCP6.0、RCP8.5 共 4 个 典型 浓 
度 路 径 情 景 (Representative concentration pathways, 
RCPs) 下 的 气候 变化 预测 结果 。 这 些 输出 结果 中 包 
括 气 温和 降水 等 气象 数据 ,时 间 跨 度 为 1850s 一 
2100 年 ,时 间 分 辨 率 有 年 、 月 、 日 3 类 。RCP6.0 和 
RCP4.5 比较 相近 ,有 旦 相当 多 的 气候 模式 没有 给 出 
RCP6.0 情 景 的 结果 ,因此 本 文 不 考虑 RCP6.0 人 情景 ， 
只 选取 RCP2.6、RCP4.5 fll RCP8.5 共 3 个 典型 浓度 
路 径 情 景 , 下 载 相应 情景 月 气温 和 降水 数据 ,共计 
39 个 ( 表 3)。 

统计 降 尺 度 方法 是 将 气候 模式 输出 的 气象 数 
据 转 换 为 流域 尺度 气象 数据 的 主要 方法 之 一 1。 
本 文 使 用 Delta 统 计 降 尺度 方法 “将 所 选取 的 各 模 
式 降 水 气温 数据 从 模式 格 网 降 尺 度 到 每 个 基站 。 
该 方法 是 直接 建立 模式 模拟 结果 与 地 面 气候 变量 
之 间 关 系 , 相 对 简单 ,并 且 随 着 GCM 模 式 结果 的 时 

空 分 辩 率 的 增加 ,其 可 靠 性 增加 。 


2 结果 与 分 析 


2.1 模型 参数 的 率 定 和 验证 

2.1.1 径流 表 4 为 3 条 河流 模型 率 定期 和 验证 期 
的 径流 评价 结果 。 率 定期 ,3 条 河流 的 多 年 平均 径 
流 模 拟 值 和 实测 值 相差 较 小 , 均 在 2% 左 右 。 验 证 
期 ,3 条 河流 的 多 年 平均 径流 模拟 值 和 实测 值 相差 
2%~5% 左 右 。 根 据 Moriasi 等 “模拟 评价 指标 , 结 
果 表 明 在 率 定 期 和 验证 期 除了 库 车 河 的 NSE 和 
PBIAS 以 外 其 余 均 达到 了 ”好 ?的 标准 。 虽 然 库 车 
河 在 率 定期 和 验证 期 的 NSE 均 为 0.63 , PBIAS 分 别 


#3 CMIPS 中 39 个 


气候 模式 基本 信息 


Tab.3 Brief overview of the 39 climate models in CMIPS 


模式 名 称 所 属国 家 4 HER 情景 
ACCESSE1-3 澳大利亚 1.25°x1.88° RCP4.5 , RCP8.5 
ACCESSEI-0 澳大利亚 1.25°x1.88° RCP4.5 .RCP8.6 
BCC-CSMI-1 中 国 2.80°x2.80° RCP2.6 , RCP4.5 , RCP8.5 
BCC-CSMI-1-m 中 国 1.12°x1.25° RCP2.6 ,RCP4.5 , RCPS.5 
BNU-ESM-3 中 国 2.80°X2.80° RCP2.6 , RCP4.5 , RCP8.5 
CanESM2 INK 2.80°x2.80° RCP2.6, RCP4.5 , RCP8.5 
CCSM4 4 美国 0.94°x1.25° RCP2.6 , RCP4.5 , RCP8.5 
CESMI-BGC 美国 0.94°x1.25° RCP4.5 , RCP8.5 
CESM1-CAM5 美国 0.94°x1.25° RCP2.6, RCP4.5 , RCP8.5 
CESMI-WACCM ”美国 1.90°x2.50° RCP2.6 ,RCP4.5 , RCPS.5 
CMCC-CESM RUM] 3.759x3.75? RCP8.5 
CMCC-CM 意大利 ”0.75°x0.75°% RCP4.5 RCP8.5 
CMCC-CMS 意大利 ”1.86°x1.875°RCP4.5 , RCP8.5 
CNRMCM 法 国 1.40°x1.40° RCP2.6 .RCP4.5 .RCP8.5 
CSIRO-Mk3-6-0 澳大利亚 1.88°x1.88° RCP2.6 .RCP4.5 .RCP8.5 
EC-EARTH ARAS 1129x125? RCP2.6,RCP4.5 .RCP8.5 
FGOALS-g2 中 国 2.81°X2.81° RCP2.6, RCP4.5 , RCP8.5 
FIO-ESM 中 国 2.792,81? RCP2.6, RCP4.5 , RCPS.5 
GFDL-CM3 美国 2.00°x2.50° RCP2.6 .RCP4.5 , RCP8.5 
GFDL-ESM2G 美国 2.00°X2.50° RCP2.6 .RCP4.5 .RCP8.5 
GFDL-ESM2M 美国 2.00°x2.50° RCP2.6 , RCP4.5 RCP8.5 
GISS-E2-H 美国 2.00°X2.50° RCP2.6  RCP4.5 , RCP8.5 
GISS-E2-H-cc 美国 2.00°x2.50° RCP4.5 
GISS-E2-R 美国 2.00°x2.50° RCP2.6 , RCP4.5 , RCP8.5 
GISS-E2-R-cc 美国 2.00°x2.50° RCP4.5 
HadGEM2-AO 韩国 1.25°x1.88° RCP2.6 , RCP4.5 , RCP8.5 
HadGEM2-CC 英国 1.25°x1.88° RCP4.5 , RCP8.5 
HadGEM2-ES 英国 1.25°x1.88° RCP2.6 , RCP4.5 , RCPS.5 
INMCM4 俄国 1.50°x2.00° RCP4.5 , RCP8.5 
IPSL-CM5A-LR 法 国 1.86°x3.75° RCP2.6 , RCP4.5 , RCPS.5 
IPSL-CM5A-MR E 1.25°x2.50° RCP2.6 ,RCP4.5 .RCP8.5 
IPSL-CM5B-LR 法 国 1.9°x3.75°  RCP4.5, RCP8.5 
MIROCS 日 本 1.40°x1.40° RCP2.6 , RCP4.5 , RCP8.5 
MIROC-ESM 本 — 2.819Q.81? RCP2.6 ,RCP4.5 , RCP8.5 
MIROC-ESM-CHEM 日 本 2.81°x2.81° RCP2.6, RCP4.5 , RCP8.5 
MPI-ESM-LR 德国 1.86°x1.88° RCP2.6 ,RCP4.5 ,RCP8.5 
MPI-ESM-MR 德国 1.86°x1.88° RCP2.6, RCP4.5 .RCP8.5 
MRI-CGCM3 日 本 1.12°x1.13° RCP2.6 ,RCP4.5 , RCPS.5 
NorESM1-M 挪威 1.90°x2.50° RCP2.6 ,RCP4.5 , RCP8.5 
Ky 6.65 和 -2.08 ,达到 模拟 效果 分 级 评价 中 “满意 ” 
标准 。 

图 2 所 示 为 模拟 的 3 条 河流 的 月 径流 与 实测 月 
径流 变化 。 从 3 条 河流 的 模拟 结果 看 ,该 模型 模拟 


的 月 径流 过 


寺 程 线 总 体 上 可 以 反应 


出 流域 实测 径流 
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表 4 SWAT-RSG 模型 在 库 车 河 \ 玛 纳 斯 河 、 库 玛 拉克 河 径 流 模拟 评价 


Tab.4 Evaluation indices for assessing runoff simulation of the SWAT-RSG in Kuqa River, Manas River and Kumaric River 


月 模拟 多 年 平均 径流 量 /10* ms 
河流 水 文 站 NSE PBIAS 率 定期 验证 期 
率 定 期 验证 期 EN 验证 期 实测 值 模拟 值 实测 值 模拟 值 
库 玛 拉克 河 协和 拉 0.66 0.85 1.55 4.08 45.50 44.80 47.30 50.00 
玛 纳 斯 河 肯 斯 瓦特 0.88 0.87 -0.65 -4.45 12.36 12.04 11.43 12.05 
库 车 河 兰 干 0.63 0.63 6.65 -2.08 2.97 3.03 3.59 3.50 
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图 2 库 玛 拉克 河 、 玛 纳 斯 河 、 库 车 河 在 模型 率 定 


期 和 验证 期 的 观测 月 径流 和 模拟 月 径流 变化 


Fig. 2 Observed and simulated monthly stream flow of Kuqa River, Manas River, 


and Kumaric River in calibration and validation periods 


的 变化 趋势 。 综 上 ,本 文 所 使 用 的 SWAT-RSG 模 型 
能 够 较 好 的 模拟 出 研究 区 径流 的 年 际 和 年 内 变化 。 
2.1.2 冰川 面积 变化 ”模拟 的 库 玛 拉克 河流 域 . 玛 
纳 斯 河流 域 和 库 车 河流 域 自 1961 年 到 2006 年 的 冰 
川 面积 分 别 退 缩 了 7.2%、19.8% 和 31.5%。2 次 冰川 
编目 “及 AsiaCryoWeb 发 布 的 TSGV1 数据 集 3 给 
出 的 冰川 面积 与 模拟 值 相差 较 小 ,分 别 退 缩 5.3%、 
19.0% 和 35.5%。 综 上 ,SWAT-RSG 模 拟 的 3 £s 
冰川 面积 变化 可 以 较 好 地 反映 各 流域 冰川 面积 和 
历史 变化 情况 。 
2.2 流域 未 来 时 期 气温 和 降水 变化 

RCP2.6, RCP4.5 和 RCP8.5 3 种 情景 下 ,各 流域 
未 来 39 个 气候 模式 集合 平均 气温 均 呈 上 升 趋势 ,日 
随 着 排放 情景 的 增强 ,气温 上 升幅 度 增 大 , 见 图 3; 


分 时 期 统计 见 表 4。3 种 情景 下 本 研究 预 估 未 来 远 
期 的 气温 变化 范围 分 别 为 1.8~1.9% 、2.5~3.0 C 、4.2~ 
6.0 ,与 张 艳 武 等 利用 CMIP5 中 29 个 气候 模式 
预 佑 的 21 世纪 末期 在 3 种 情景 下 中 国平 均 气 温 分 
别 升 高 2.0% 3.3 *C [16.3 % 相 比 偏 低 ,其 主要 原因 
iach a ionic 高 海拔 地 区 。 
种 情景 下 , 随 着 排放 情景 升 高 , 预 估 的 各 流域 
ER e 图 3) ,但 分 
时 期 统计 各 流域 未 来 时 期 降水 较 历 史 时 期 均 有 增 


加 ,增加 幅度 较 小 ( 表 5)。 从 表 5 可 以 看 出 ,降水 变 
化 在 不 同情 景 下 及 流域 间 有 些许 差异 。 


2.3 流域 未 来 径流 变化 
2.3.1 年 际 变化 TERCP2.6, RCPA.5 和 RCP8.5 人 情 
景 下 ,未 来 时 期 3 个 流域 的 径流 量 较 历 史 时 期 均 有 
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注 : 实 线 为 39 个 气候 模式 集合 平均 值 ;阴影 部 分 为 39 个 气候 模式 集合 的 最 大 值 .最 小 值 表示 的 不 确定 性 区 间 。 下 同 。 


图 3 RCP2.6.RCP4.5.RCP8.5 tit P 1961 一 2100 年 库 车 河 、 玛 纳 斯 河 、 库 玛 拉克 河 气 温和 降水 年 际 变 化 趋势 
Fig. 3 Inter-annual change trends of temperature and precipitation in Kuqa River, Manas River, 
and Kumaric River during 1961—2100 under RCP2.6, RCP4.5 and RCP8.5 scenarios 


增加 (图 4、 表 5), 这 与 Lutz 等 % 预 估 的 未 来 亚洲 高 
海拔 地 区 5 个 流域 上 游 径流 变化 趋势 一 致 。 冰 融 水 
变化 在 不 同时 期 ,不同 情景 及 流域 间 存 在 差异 , 库 
车 河和 库 玛 拉克 河 在 未 来 远 期 RCP8.5 情境 下 变化 
最 大 为 -65.1% 和 208.2% ; 玛 纳 斯 河 则 在 未 来 远 
RCP2.6 情 境 下 变化 最 大 ,减少 了 37.9%( 表 5), 受 冰 
融 水 影响 各 流域 未 来 2 个 时 期 的 径流 量 较 历史 时 期 
的 变化 也 存在 差异 ,对 于 冰川 覆盖 率 较 大 , 冰 融 水 
贡献 大 的 库 玛 拉克 河流 域 , 总 体 上 呈 先 增加 后 减少 
的 趋势 ,这 与 一 些 人 研究 学 者 ”” 给 出 的 冰川 覆盖 比 
例 越 高 的 流域 ,其 未 来 径流 可 能 增加 得 越 多 的 结论 
一 致 。 玛 纳 斯 河流 域 的 径流 在 3 种 情景 下 均 呈 先 增 
加 后 减 小 的 趋势 , 且 径流 的 变化 趋势 与 降水 的 增多 
趋势 不 一 致 。 而 对 于 冰川 覆盖 率 极 小 , 主要 由 降水 
补给 为 主 的 库 车 河流 域 , 在 3 种 情景 下 径流 均 有 增 
加 。 综 上 ,在 冰川 覆盖 率 小 的 库 车 河流 域 ,径流 与 
降水 变化 趋势 一 致 ,而 在 冰川 覆盖 率 较 大 的 玛 纳 斯 
河和 库 玛 拉克 河 径 流 呈 先 增加 后 减少 或 持续 增加 


趋势 ,与 降水 变化 趋势 不 完全 一 致 。 
2.3.2 年 内 分 布 在 3 种 情景 下 ,降水 变化 差异 小 ， 
但 在 RCP8.5 情景 下 气温 升 高 幅度 最 大 ,引起 的 径 
流 和 雪 冰 融 水 变化 也 更 为 显著 。 因 此 ,本文 着 重 
探讨 在 RCP8.5 情景 下 不 同 冰川 覆盖 率 的 流域 雪 冰 
融 水 的 季节 变化 ,从 而 揭示 径流 季节 变化 的 主要 
预 估 各 流域 未 来 春季 径流 量 较 历史 时 期 均 有 
增加 ,这 种 变化 主要 与 春季 融雪 峰值 时 间 提 前 、 融 
雪 量 增 大 有 关 。 夏 季 各 河流 的 径流 量 峰值 时 间 无 
明显 偏 移 ,但 是 径流 峰值 量 在 不 同 冰川 覆盖 率 的 流 
域 发 生 不 同 程 度 的 变化 (图 5)。 冰 川 覆 羡 率 最 大 的 
库 玛 拉 克 河 夏季 融雪 量 在 未 来 2 个 时 期 均 有 减少 ， 
但 是 其 最 大 径流 量 受到 冰 融 水 在 夏季 显著 增加 的 
影响 ,径流 量 较 历 史 时 期 均 有 增加 。 值 得 注意 的 是 
该 流域 近期 的 最 大 径流 量 大 于 远 期 ,与 冰 融 水 变化 
不 一 致 ,推测 其 原因 在 于 该 情景 下 气温 增加 明显 导 
致 该 流域 夏季 蒸 散 加 大 。 
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表 5 研究 区 未 来 近期 和 远 期 39 个 气候 模式 集合 的 气温 、 降 水 、 径 流 、 冰 融 水 均值 较 历 史 时 期 变化 


Tab.5 Changes of temperature, precipitation, streamflow, and ice melt on averages of the 39 climate models for the two 


future periods relative to the historical period in the study area 


Him 5E 库 车 河 玛 纳 斯 河 库 玛 拉克 河 
气温 /C ”降水 /% 冰 融 水 /9% 径流 /% ”气温 /C KEk% 冰 融 水 /% 径流 /% ”气温 /C ”降水 /% 冰 融 水 /% 径流 /% 
平均 值 1.7 4.7 -5.7 6.4 1.5 3.7 25.0 10.1 1.7 33 111.8 28.3 
RCP2.6 WAW 2.9 23.6 312 51.6 2.9 18.1 161.5 41.8 2.6 155 212.6 56.2 
(近期 ) ”最 小 值 ^ 06  -104 -40.7  -257 0.0 -97  -430 -88 06 10.5 4.6 82 
标准 差 0.6 74 15.1 15.0 0.7 5.9 61.9 13.9 0.5 6.1 523 13.4 
平均 值 1.8 3.9  -646 7.8 1.7 47  -319 2.9 1.9 3.8 40.5 12.9 
RCP26 ”最 大 值 3.6 280 -49.1 72.8 3.1 24.9 12.1 32.8 3.5 17.2 111.5 35.9 
( 远 期 ”最 小 值 -03  -113 -77.1  -231  -02 -132  -516  -181 -02 -9.1 -250  -111 
标准 差 0.8 9.4 73 20.5 0.9 74 13.6 10.7 0.8 7.0 33.5 11.4 
平均 值 1.6 4.4 -2.1 7.4 1.4 3.8 23.5 9.8 1.7 3.2 1211 30.5 
RCP4.5 ”最 大 值 27 21.9 26.4 43.2 2.4 17.7 156.5 39.7 2.7 129 226.8 62.9 
GEW) ”最 小 值 —-13 -9.5  -279 -212  -01  -137 -49.0 -7.2 -13 -12.0 10.8 6.5 
标准 差 0.7 7.0 12.1 14.0 0.6 6.3 58.6 12.6 0.7 5.8 52.2 11.7 
平均 值 3.0 63 | -600 12.4 2.5 58  -220 6.6 3.0 44 — 1034 27.0 
RCP4.5 RAE 4.9 29.2 -33.3 73.8 4.8 25.4 13.9 40.8 4.9 18.1 154.2 48.0 
( 远 期 ) ”最 小 值 Ll -128 -803  -268 -02  -130  -608  -129 Ll — -12.8 24.0 6.4 
标准 差 0.9 8.6 11.4 19.4 II 7.6 20.4 11.8 0.9 74 31.8 10.1 
平均 值 2.0 4.2 3.4 5.8 1.6 3.7 34.5 11.4 2.0 31 1429 34.8 
RCP&5 ”最 大 值 3.1 27.8 33.5 60.3 3.1 216 178.3 46.0 3.0 183 251.4 67.9 
(近期 ) ”最 小 值 07 -176 -213 -239  -04  -123 -47.1 -7.1 0.8 -19.7 33.5 6.8 
标准 差 0.5 7.8 12.4 15.1 0.8 6.2 71.3 14.9 0.5 7.0 54.8 13.0 
平均 值 5.2 6.0 -65.1 14.3 4.2 5.9 -0.3 9.8 52 25 2082 47.1 
RCP8.5 RAW 7.5 405  -282 120.4 73 39.1 77.8 64.7 7.4 273 260.0 78.0 
( 远 期 ) ”最 小 值 26  -247 -883 -373 -01 -194 -492  -174 25  -253 1034 25.8 
标准 差 1.2 12.6 15.0 29.1 1.9 10.4 37.8 16.4 1.2 10.9 32.5 11.8 


24 预 估 结 果 的 不 确定 性 

流域 水 文 模拟 的 不 确定 性 主要 来 源 于 输入 数 
据 、 水 文 模型 结构 和 水 文 模型 参数 3 个 方面 。 本 
研究 已 经 利用 历史 观测 水 文 资料 和 冰川 变化 数据 
对 水 文 模型 的 模拟 能 力 进 行 了 相关 评价 ,达到 了 满 
意 的 标准 。 而 气候 模式 输出 的 未 来 气象 数据 作为 
水 文 模 型 的 驱动 因子 ,是 预 佑 流域 水 文 过 程 结果 不 
确定 性 的 重要 来 源 “”。 前 人 的 研究 表明 多 个 气候 
模式 集成 均值 和 单个 模式 相 比 ,能 更 好 的 反映 出 观 
测 数据 的 变化 趋势 ”。 

本 文集 中 39 个 气候 模式 ,在 年 和 月 尺度 上 进行 
模拟 预 估 ,由 不 同 气 候 模 式 给 出 的 气温 和 降水 输入 
水 文 模型 引起 流域 模拟 结果 变化 的 全 部 变化 范 
围 。 在 年 尺度 上 (图 4) ,冰川 覆盖 率 越 大 的 流域 , 预 
佑 径流 不 确定 性 变化 区 间 受 冰 融 水 影响 越 大 。 而 


注 :平均 值 . 最 大 值 和 最 小 值 对 应 的 正 值 表示 增加 , 负 值 表示 减少 ; 表 


FP 的 气温 为 绝对 变化 ,降水 径流 和 冰 融 水 均 为 相对 变化 。 


冰 融 水 的 不 确定 性 区 间 总 体 表现 为 未 来 近期 大 于 
未 来 远 期 ,在 未 来 近期 RCP2.6 和 RCP8.5 这 2 个 极端 
的 情景 下 的 不 确定 区 间 大 于 RCP4.5 情景 , 旦 冰川 覆 
盖 率 较 大 的 流域 冰 融 水 不 确定 性 区 间 较 大 。 图 6 给 
出 了 39 个 气候 模式 输出 的 气象 数据 驱动 水 文 模型 
预 估 流 域 未 来 近期 和 远 期 较 历 史 时 期 年 内 径流 变 
化 的 不 确定 性 区 间 。 预 估 的 各 流域 夏季 径流 变化 
不 确定 性 区 间 明 显 大 于 其 它 季 节 及 年 平均 值 (图 6)， 
这 主要 与 夏季 冰 融 水 发 生 明 显 变化 有 关 ( 图 $)。 且 
不 同 冰川 覆盖 率 的 流域 其 不 确定 性 区 间 有 明显 差 
异 , 冰 川 覆 盖 率 越 大 , 预 估 的 夏季 径流 变化 的 不 确 
定性 区 间 在 未 来 2 个 时 期 和 不 同情 景 下 差异 越 大 
(图 6)。 且 夏季 39 个 气候 模式 预 估 径流 变化 的 离散 
程度 增加 ,可 能 导致 极端 水 文 事件 也 增加 ,尤其 是 
在 RCP8.5 情景 下 各 流域 径流 变化 不 确定 性 区 间 大 
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图 4 RCP2.6.RCP4.5,RCP8.5 tit F 1961 一 2100 年 库 车 河 ` 玛 纳 斯 河 库 玛 拉克 河 径流 和 冰 融 水 年 际 变化 趋势 
Fig.4 Inter-annual change trends of streamflow and ice melt in Kuqa River, Manas River, 
and Kumaric River during 1961—2100 under RCP2.6, RCP4.5 and RCP8.5 scenarios 


一 径流 1966—19954E — 雪 融 水 _1966 一 1995 年 ”一 冰 融 水 _1966 一 1995 年 
-- 径流 2016—20454E  -- 雪 融 水 2016 一 2045 年 ” -- 冰 融 水 _2016 一 2045 年 
- 径流 2066 一 2095 年 ”… 雪 融 水 _2066 一 2095 年 “” … 冰 融 水 _2066 一 2095 年 
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注 : 正 误差 棒 代 表 2016 一 2045 年 39 模 式 标准 差 , 负 误 差 棒 代 表 2066 一 2095 年 39 模式 标准 差 。 

图 5 RCP8.5 情 景 下 库 玛 拉克 河 , 玛 纳 斯 河 , 4 4-30 1966—1995 ^F. 2016—2045 年 .2066 一 2095 年 雪 冰 产 流 及 径流 年 内 分 布 
Fig.5 Intra-annual distribution of snow melt, ice melt and stream flow in Kuqa River, Manas River and 
Kumaric River for the three periods of 1966—1995s, 2016—2045s, 2066—2095s under RCP8.5 scenario 


于 另外 2 个 情景 下 的 径流 变化 不 确定 性 区 间 。 RCP8.5 情景 下 , 预 估 不 同 冰 川 履 盖 率 的 流域 未 来 雪 
冰 融 水 及 径流 的 变化 ,讨论 径流 对 未 来 气候 变化 的 
3 结论 


响应 及 不 确定 性 ,得 出 以 下 结论 : 
本 文 基于 SWAT-RSG 水 文 模型 ,将 气候 模式 数 (1) 随 着 未 来 气温 和 降水 增加 ,各 流域 的 径流 
据 降 尺 度 处 理 ,驱动 水 文 模型 ,在 RCP2.6、RCP4.5 和 — 均 呈 增加 趋势 ,但 是 不 同 冰川 覆盖 率 的 流域 其 径流 
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红 等 : 多 模式 预测 气候 变化 及 其 对 雪 冰 流域 径流 的 影响 
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证:0.95 表 示 预 估 的 径流 在 95% 分 位 上 的 变化 值 ;0.05 表 示 预 估 的 径流 在 5% 分 位 上 的 变化 值 。 


图 6 RCP2.6、RCP4.5 .RCP8.5 情 景 下 39 个 气候 模式 驱动 模型 预 估 流 域 2016 一 2045 年 和 
2066 一 2095 年 较 1996 一 1995 年 径流 变化 箱 线 图 
Fig.6 Boxes of runoff changes in the basin during 2016—2045 and 2066—2095 compared with 1996— 1995 
by the 39 climate model-driven models under RCP2.6, RCP4.5 and RCP8.5 scenarios 


的 年 内 和 年 际 变化 在 未 来 时 期 会 发 生 较 大 的 差 
别 。 尤 其 是 在 冰川 覆盖 率 大 的 库 玛 拉克 河流 域 径 
流 主要 受到 雪 冰 融 水 变化 影响 ,对 气候 变化 响应 更 
显著 , 随 着 排放 情景 增加 径流 量 显著 增加 ,在 年 内 
夏季 径流 量 较 历 史 时 期 也 显著 增加 。 冰 川 覆 盖 率 
很 小 的 库 车 河流 域 ,径流 与 降水 变化 趋势 一 致 ,在 
年 内 该 河流 夏季 径流 量 则 表现 出 减 小 的 变化 。 

(2) 39 个 气候 模式 集合 驱动 水 文 模型 ,集合 平 
均值 反映 了 绝 大 部 分 气候 模式 数据 可 能 导致 的 未 
来 水 文 过 程 的 变化 ,波动 区 间 则 反映 了 该 区 域 水 文 
过 程 潜在 的 不 确定 性 。 冰 川 覆 盖 率 越 大 的 流域 其 
径流 预 估 的 不 确定 性 受 冰 融 水 影响 越 大 , 且 在 
RCP2.6 和 RCP8.5 极 端 排放 情景 下 不 确定 性 区 间 大 
于 RCP4.5 情 景 。 在 冰川 覆盖 率 大 的 流域 其 径流 不 


确定 区 间 变 化 与 冰 融 水 变化 一 致 , 均 表 现 为 未 来 近 
期 的 不 确定 区 间 大 于 未 来 远 期 。 对 于 冰川 覆盖 率 
很 小 的 流域 ,径流 变化 的 不 确定 性 区 间 在 未 来 2 个 
时 期 逐渐 增 大 ,与 该 流域 降水 不 确定 性 变化 有 关 。 

预 佑 结果 为 39 个 气候 模式 集合 ,包含 了 绝 大 部 
分 气候 模式 ,减少 了 由 单个 或 几 个 气候 模式 预 佑 的 
气温 和 降水 数据 驱动 水 文 模型 带 来 的 不 确定 性 ,但 
由 于 气候 模式 本 身 的 不 确定 性 ,以 及 不 同 气 候 模 式 
选取 集合 方式 对 预 佑 结果 也 会 产生 不 同 程度 的 影 
Ul] ,应 做 进一步 不 确定 性 分 析 人 研究 ,提高 未 来 气候 
变化 对 流域 水 文 过 程 影响 预 佑 中 的 准确 性 。 
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Predicting climate change and its impact on runoff in 
snow-ice basin with multi-climate models 
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Abstract: The snow-ice melt water in glacier-covered basins plays an essential role in regulating runoff. The im- 
pact of climate change on the snow-ice melting process and quantity change, river runoff process and quantity 
change, which is related to the glacier area in the basin. Using climate models’ meteorological data of CMIPS5 
driven watershed hydrological model, simulated runoff response to climate change in Tianshan areas (Kumaric 
River, Manas River, and Kuga River, Xinjiang, China). The results show that with the temperature and precipita- 
tion increasing in the future, three river basins’ snow melt water is both increasing, but ice melt water changes af- 
fected by glacier coverage. The change of runoff was mainly affected by the increase in precipitation and the 
change in snow and ice melt water; the runoff increased by 5.8%-14.3% (Kuga River), 2.9%-11.4% (Manas River), 
and 12.9%-47.1% (Kumaric River) under three representative concentration pathways (RCPs). It was estimated 
that the variation range of runoff uncertainty was affected by ice melting. From the projected intra-annual distribu- 
tion of runoff, the snow and ice melt of the three basins, the spring runoff of the river basins increased compared 
with the historical period due to the snow melt time ahead and snow melt volume. In the summer, the peak runoff 
of the Kuqa River and Manas River decreased due to ice melt water, whereas the peak value of runoff in Kumaric 
River increased. In addition, the estimated uncertainty range of summer runoff in the study basins in the summer 
season is significantly larger than those of other seasons. 
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